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Abstract 



We study the synchrotron-like gluon emission in AA-collisions from fast partons due to 
interaction with the coherent glasma color fields. Our results show that for RHIC and 
LHC conditions the contribution of this mechanism to parton energy loss is much smaller 
than the radiative energy loss in the plasma phase. 

jjn ! 1 Introduction 

It is widely believed that the observed large suppression of high-j»T hadrons in AA- 
collisions at RHIC and LHC (the so-called jet quenching (JQ)) is due to parton energy 
loss in the quark-gluon plasma (QGP). The observation of this phenomenon and success- 



ful hydrodynamic description of hadron production at low are the strongest arguments 



for formation of the QGP in heavy-ion collisions at RHIC and LHC. In pQCD partons 
lose energy mostly through the induced gluon emission due to multiple scattering in the 
(N) ■ QGP [1, 2, 3, 4, 5, 6], and a small fraction of the energy loss comes from the elastic 
collisions [7, 8]. At RHIC and LHC, the entropy of the QGP required to describe the 
observed suppression of high-p-r hadrons, within the pQCD calculations of the radiative 
and collisional energy losses [9, 10, 11], turns out to be qualitatively consistent with that 
obtained within the hydrodynamic simulations from the low-p^ data. This consistency 
looks very encouraging. However, one must keep in mind that the theoretical uncertain- 
ties both in the pQCD calculations of the JQ and in the hydrodynamic simulations are 
rather large. On the JQ side there remain open questions on additional mechanisms of 
the parton energy loss which can modify the picture with the dominating radiative energy 
loss in the QGP phase. 

One of the potentially important energy loss mechanisms is the synchrotron-like gluon 
emission in the strong classical color fields in the initial pre-equilibrium stage of AA- 
collisions [12] predicted within the Color Glass Condensate (CGC) model [13, 14, 15, 
16]. This stage (termed glasma [17]) is now under active investigation (see [18, 19] and 
references therein). However, its role in parton energy loss still has not been studied. 
The evaluation of the classical Yang-Mills fields in AA-collisions within the CGC model 
shows that just after the collision of the Lorentz contracted nuclei a system of the color 
flux tubes with the longitudinal boost invariant color electric and color magnetic fields 
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(with \E Z \ \B Z \) should be produced [20]. The typical transverse coherence length for 
the color fields in this phase is ~ 1/<5 S , where Q s (~ 1 — 1.5 GeV for RHIC and LHC 
conditions [21]) is the saturation scale of the nuclear parton distributions. At the proper 
time r ~ l/Q s the strength of the glasma electric and magnetic longitudinal fields is 
gE z ~ gB z ~ Q 2 S [20]. The transverse fields, which are absent at r = 0, are generated 
at later times and become close to the longitudinal ones at r > 1/Q S [20, 22]. At such 
times the glasma energy density decreases oc 1/r. At later times the thermalization of 
the glasma color fields should lead to formation of the equilibrium QGP, but the detailed 
mechanism for the thermalization remains unclear. Qualitative analyses say that the 
glasma thermalization goes probably via instabilities of the color flux tubes rising quickly 
at T ^^/Qs [23, 22, 24]. They should lead to a fast randomization of the boost invariant 
glasma color fields at the time scale about 2-3 units of 1/Q S . This thermalization time 
agrees qualitatively with the hydrodynamic fits of the data on AA-collisions which favor 
the thermalization at r ~ 0.4 — 1 fm [25, 26]. 

For RHIC-LHC conditions the typical Lorentz force acting on a fast parton crossing the 
glasma turns out to be very large ~ Q 2 S ~ 5 — 10 GeV/fm. It is about 10-20 times that for 
the Debye screened color center in the QGP ~ a s m 2 D (if one takes a s ~ 0.3 and m D ~ 0.5 
GeV). It raises a natural question about the impact of the glasma on the radiative energy 
loss. If one ignores the interaction of fast partons with the pre-equilibrium phase at all, 
then for jets with energies < 100 GeV the parton showering at r < t ~ 0.5 — 1 fm can be 
viewed as almost pure DGLAP shower, and the interference of the DGLAP stage with the 
induced gluon emission in the bulk of the QGP at r > t should be relatively small [10]. 
Thus in this energy range the question of the glasma role in the JQ is equivalent to the 
question to what extent the glasma modifies the DGLAP gluon emission. The purpose of 
the present paper is to give a qualitative estimate of this effect. 

2 Formulation of the model and basic formulas for 
the synchrotron-like gluon emission 

An accurate calculation of the synchrotron-like gluon emission due to interaction of fast 
partons with the glasma is presently impossible. It would require a detailed understanding 
of the decay and thermalization of the glasma color flux tubes which remain unclear. 
However, since the strength of the color fields in the flux tubes decreases quickly, the 
transverse momentum which a fast parton gets in the glasma should mostly come from 
interaction with the coherent color field of the first crossed color flux tube. The effect of 
the random transverse momentum kicks which a fast parton undergoes in the glasma at 
later times should be small due to weakness of the fields and the destructive interference of 
the chaotic contributions from different color tubes. For the above reasons as a first step 
in understanding the influence of the glasma on parton energy loss it seems reasonable to 
consider a simple model with a uniform time-dependent color field which acts only for a 
limited range of r about 2-3 units of l/Q s . 

So we will consider the synchrotron-like gluon emission from a fast parton (we choose 
the z-axis along the initial parton momentum) in a slab of thickness L ~ (2 — 3) /Q s with 
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the transverse chromoelectric, E , and chromomagnetic, H a , fields (hereafter a denote the 
color index). We will ignore the longitudinal (along the jet) components of the color fields 
since similarly to the photon emission in QED their effect should be small for relativistic 
partons. For simplicity we assume that both the electric and magnetic fields are oriented 
to the same direction in the color space. It is enough to consider the color fields with 
nonzero components only in the Cartan subalgebra, i.e., for a = 3 and a = 8 for the 
SU (3) color group. 

We present our formulas for a fast quark produced in a hard process at central rapidity 
y = (so the initial quark momentum and the z-axis of our coordinate system are 
perpendicular to the AA-collision axis). We take z = for the production point. The 
starting point of our analysis is similar to the case of the synchrotron-like energy loss 
spectrum for a fast parton propagating in an infinite uniform color field addressed in [27]. 
We write the S'-matrix element of the q — > gq' transition as (we omit the color factors 
and indices) 

(gq'\S\q) = -ig J d 4 y^(y) 7 M;(y)^(y) , (I) 

where ip qiq > are the wave functions of the initial quark and final quark, is the wave 
function of the emitted gluon. 

For our choice of the external color field (which has only a = 3 and a = 8 nonzero 
components) it does not change the quark color state. However, it is not the case for 
gluons, which are rotated in the color space as they move in the color field. The interaction 
of the gluons with the background color field can be diagonalized by introducing the gluon 
fields having definite color isospin, Qa, and color hypercharge, Qb, (we will describe 
the color charge by the two-dimensional vector Q = (Qa,Qb))- In terms of the usual 
gluon vector potential, A a , the diagonal color gluon states read (the Lorentz indices are 
omitted) X — (Ai + iA 2 )/y/2 (Q = (-1, 0)), Y = (At + iA 5 )/V2 (Q = (-1/2, -y/Z/2)), 
Z = (A 6 + iA 7 )/y/2 (Q = (1/2, — v/3/2)). The neutral gluons A = A 3 and B = A 3 with 
Q = (0, 0), to leading order in the coupling constant, do not interact with the background 
field at all, and the emission of these gluons are similar to the photon radiation in QED. 

We evaluate the matrix element (1) within the small angle approximation, i.e., we as- 
sume that for each parton its transverse momentum is small compared to the longitudinal 
one. We write each quark wave function in the form 

^H e *P'-^-*>W .P). (2) 

where the four-vector y is defined as y^ = (t, p, z) (hereafter the bold vectors denote the 
transverse vectors), A is the quark helicity, u\ is the Dirac spinor operator. In the small 
angle approximation the z-dependence of the transverse wave functions 0« is governed by 
the two-dimensional Schrodinger equation 

* -q z = { 2E_ q - + U i (z,p)\(j) i (z,p) (3) 

with the potential 

U i (z,p) = gQ i n [G n (z,p)-G 3 n (z,p)}. (4) 
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Here G£ denotes the external vector potential (the superscripts are the Lorentz indexes 
and n = 1, 2 correspond to the A and B color components in the Cartan subalgebra), Q l n 
is the quark color charge. 

We take the external vector potential in the form G° = — p • E n , G n = 0, and 
G\ = [H„ xpf, where E n and H n are the electric and magnetic fields. In this case the 
kinetic term in (3) does not contain the external vector potential at all, and the potential 
Ui can be written as 

U i (z,p) = -F i -p, (5) 

where Fj is the Lorentz force. The wave function of the emitted gluon can be represented 
in a similar way. In our numerical computations for the quark mass we take m q = 0.3 
GeV. However, the value of the quark mass is not very important. For the gluon mass 
we take m g = 0.75 GeV. This value was obtained from the analysis of the low-x proton 
structure function F 2 within the dipole BFKL equation [28]. It agrees well with the 
natural infrared cutoff for gluon emission m g ~ l/R c , where R c ~ 0.27 fm is the gluon 
correlation radius in the QCD vacuum [29]. 
The solution of (3) can be taken in the form 

<i>i(z,p)=exp{ip i (z)p-^J o dz'[p 2 {z')+m 2 q ]^. (6) 

Here the transverse momentum Pi(z) is the solution to the classical parton equation of 
motion in the impact parameter plane 



By using (1), (2), (6) one can obtain 



(gq'\S\q) = -ig(2Ttf5(uj + Eg, - E q ) / dzV(q(z), {\})5(p g (z) + p q ,(z) - p q {z)) 

J 




p 2 g (z')+m 2 g i pl(z') + m 2 q p 2 Q (z')+m 2 Q ]} ^ 



_ _JL + --.--. - — 
loo ^Eqi 2E q 



where u is the gluon energy, V is the vertex factor, {A} is the set of the parton helicities. 
For the transition conserving quark helicity V(q, A) = — iC[2X q x + (2 — x)X g ][q x (z) — 



i\ g q y (z)]/[xy2(l — x)], and for the spin-flip case V = ixm q C(2X q X g + 1)/ y2(l — x). Here 
q(z) = p g (z)(l —x) — p q r(z)x, x = oj/E q , and C = X^Xa/Z is the color factor (i, f are the 
color indexes of the initial and final quarks, Xa is the color wave function of the emitted 
gluon). Due to the color charge conservation F q = F g + F q/ . For this reason, the argument 
of the momentum 5-function on the right-hand part of (8) does not depend on z, and can 
be replaced by p+ + pj — p^j", where pf = pj(oo). Then, noting that the expression in the 
square brackets in (8) can be rewritten as [q 2 (V) + e 2 ]/2M with e 2 = m 2 x 2 + (1 — x)m 2 
and M = E q x(l — x), the matrix element (8) can be rewritten as 

(gq'\S\q) = -*(2tt) 3 5( W + E q , - E q )6(p+ + pj - p+)T , (9) 
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g J dzV(q(z), {A}) exp < i J dz' 



qV) + e 2 
2M 



(10) 



With the help of the standard Fermi golden rule the gluon distribution in terms of the 
amplitude T can be written as 



\T\ 



dN 



dudq 8(27r) 3 £ 3 :r(l - x) ' 



Hereafter we do not show the averaging over the initial and summing over the final 
color/spin states. In the absence of the external color field the momentum q(z) in (10) 
does not vary with z and equals to the final momentum at z = oo (we will denote it q). 
In this case (10) gives the vacuum amplitude 



T v = gV(q,{\}) 



2iM 
q 2 + e 2 



(12) 



Substituting (12) into (11) we obtain the usual vacuum LO pQCD distribution for the 
q — )■ gq splitting (hereafter we neglect the spin-flip part of the vertex (oc m q ) which gives 
a negligible contribution) 



dN„ 



2C F a s 



dudq n 2 xE q V 2 / (q 2 + e 2 ) 2 ' 

For a nonzero external field we write the amplitude T as a sum 

T = T V + T S , 



(13) 



(14) 



where T s describes the correction due to the external field. Since the momentum q varies 
with z only at z < L, from (10) one can obtain 



T s = g [ L dzV{<i(z),{\})exp 
Jo 



% I dz' 



q 2 (^)+6 2 
2M 



- (q(^) q) 



(15) 



In terms of T VjS the synchrotron-like correction to the LO gluon spectrum (13) reads 



dN s 2Re(T v T s *) + \T S 



dudq 8(27r) 3 E3a;(l - x) 



(16) 



The glasma correction to the cu-spectrum can be obtained from (16) by integrating over 
the transverse momentum q. Note that contrary to the LO vacuum spectrum (13), which 
oc 1/q 2 at large q 2 and without a constraint on q 2 gives a logarithmically divergent 
cu-spectrum, the q-integral for the correction term is convergent at large q 2 . 

The w-spectrum can also be obtained within the light-cone path integral (LCPI) ap- 
proach [2] formulated in the impact parameter space. The LCPI formalism has originally 
been developed for gluon emission induced by multiple parton scattering. But it applies to 
the synchrotron-like gluon emission as well. In this case the LCPI gives for the cu-spectrum 



dN„ 



du 



/•oo roo 

2Re / dzi dz 2 g[!C(p 2 ,Z2\p 1 ,z 1 ) - lC (p 2 ,z 2 \p 1 ,z 1 )] 

JO Jz! 



Pl = P 2 =0 



(17) 
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Here g is the vertex operator given by 



+ a A. 9 (18) 



£3x 3 (l -x) 2 V 2 ; dp 1 dp 2 ' 

/C is the Green's function of the Schrodinger equation with the Hamiltonian 



"--wA^j-'-^Wr (19 > 

and /Co is the Green's function for the Hamiltonian (19) with f = 0. The Green's function 
for the Hamiltonian (19) is known explicitly (see, for example, [30]) 

M 

fc(P 2 > 32|Pi, zi) = 2niAz ex P [iSd] , (20) 
where Az = z 2 — z 1 and S c i is the classical action for the Hamiltonian (19) given by 

Aze 2 M \, 2 f z * 

^cl — WW + 



2M 2Az 

+ if Z dtpi ' m{z2 ~ t] ~^ Z dt L dsm ' i{s){z2 ~ t){s ' zi) 



(P 2 - Pi? + ^ /J' dtp 2 ■ f(t)(t - Zl ) 



(21) 



One can show that the LCPI formula (17) being rewritten in the momentum space repro- 
duces the w-spectrum corresponding to (16). Note that the subtraction of the /Co term in 
(17) corresponds to subtraction of the vacuum part from (11) for getting the synchrotron 
correction (16). 

Note however that for RHIC and LHC conditions in our approach the energy loss 
spectrum can not be calculated accurately at gluon energy uj < 5 GeV. Our numerical 
calculations show that for such soft gluons the large angle emission becomes important and 
the integration over the transverse momentum can not be performed accurately using the 
formulas based on the small angle approximation. For this reason our results on the energy 
loss can only be treated as qualitative estimates, rather than quantitative predictions. 

We presented the formulas for the q — > gq' transition. For the purely gluonic g — >■ gg 
transition the calculations are similar. 



3 Numerical results 

For numerical calculations we should fix the ^-dependence of the Lorentz force. As we 
have said, we assume that the external field acts a finite time/length ~ 1/Q S - In our 
numerical calculations we take L = 2/Q s . We take Q s — 1 for RHIC and Q s = 1.4 for 
LHC. This gives L(RHIC) « 0.4 fm and L(LHC) w 0.28 fm. To fix the z-dependence 
of the Lorentz force, which a fast parton feels at r = z < L, we use the r-dependence 
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of the glasma energy density obtained in the lattice simulations by Lappi [20]. For our 
calculations we need only the field components transverse to the initial parton momentum. 
At t <C l/Qs, when the electric and magnetic fields are almost parallel to the AA-collision 
axis [18, 19], the Lorentz force acting on a fast parton is purely transverse to the parton 
momentum. At such times for a unit color charge F 2 = 2g 2 e, where e = (E 2 + H 2 )/2 
is the color field energy density. However, this relation is invalid at r > 1/Q S when the 
contribution of the transverse (to the AA-collision axis) components of the color fields to 
the energy density becomes approximately equal to that from the components along the 
beam axis [20] . Since only half of these transverse fields squared contribute to the Lorentz 
force squared which we need (transverse to the jet direction) we can write in this regime 
F 2 = g 2 3e/2. We will use this relation to determine the Lorentz force from the energy 
density obtained in [20]. 

The results of [20] were presented in terms of the coupling constant g and the mass 
parameter /i of the CGC model [13, 14] related to the saturation scale by the relation 
Q s Ki Qg 2 fi/4ir. Extrapolating to the continuum limit Lappi obtained e(t = l/g 2 [i) = 
0.2Q(g 2 fx) 4: /g 2 , and at r > \jg 2 \x the energy density behaves as ~ 1/r. We use the above 
value of e{t = l/g 2 fi) to fix the normalization, and the r-dependence of e was obtained 
using the results for g 2 e/(g 2 fi) A presented in Fig. 3 of [20]. In [20] the glasma energy 
density has been calculated taking g = 2, and g 2 fi = 2 GeV for RHIC and g 2 jj, = 3 GeV 
for LHC. Note that it gives the values of the glasma energy density at r = l/g 2 /i which 
match very well to the QGP energy density extracted within the Bjorken model with the 
longitudinal expansion [31] using the entropy/multiplicity ratio dS / dy j dN ch / dr] = 7.67 
obtained in [32] and experimental multiplicities for the central Au + Au collisions at 
v/s = 200 GeV at RHIC and Pb + Pb collisions at yfs = 2.76 TeV at LHC. In terms of 
the QGP temperature it corresponds to T « 300 MeV for RHIC, and T « 400 MeV for 
LHC at r = 0.5 fm. We use these parameters to fix the Lorentz force. For the coupling 
constant in the amplitude of the q — >■ qg transition (10) we also take g — 2, i.e. a s « 0.318. 
Our numerical evaluations show that the results for the external fields in the A and B 
color states are very similar. Below we present the gluon emission spectra defined as a 
half-sum of the A and B contributions. 

To illustrate the transverse momentum dependence of the gluon emission in Fig. 1 we 
plotted the distribution dN s /dudq scaled by its value at q = for q along the x and y 
axes (we take the x-axis parallel to the Lorentz force, and the y-axis perpendicular to the 
Lorentz force) for a quark with energy E = 50 GeV and u — 2, 5, 10, and 25 GeV for 
RHIC and LHC conditions. The (^-dependence in Fig. 1 is shown for the distribution 
averaged over two directions of the Lorentz force. Fig. 1 shows that the glasma correction 
has a complicated q-dependence with a region where it is negative. The appearance of 
the negative dN s / dud(\ shows that a significant contribution comes from the interference 
of T v and T s amplitudes in (16). This says that the finite-size effects are very important 
for RHIC and LHC conditions (see below). From Fig. 1 one sees that the q-distribution is 
rather broad, and large angle gluon emission is clearly important at uj < 5 GeV. However, 
our approach, based on the small angle approximation, becomes inapplicable at the angles 
> 1. For this reason for the ^-integrated energy loss spectrum dN s /du our approach 
can give only a rough estimate of the magnitude of the glasma contribution at u> < 5 
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GeV. The observed significant contribution of the large angle gluon emission is somewhat 
unexpected. Indeed, one might naively expect that the typical transverse momenta should 
be ~ Q s (since the Lorentz force is ~ Q 2 S and the path length is ~ 1/Q S ). 

We evaluated the dN s /du for two different prescriptions: (a) without any constraint 
on the transverse momentum (ignoring the fact that our formulas are inapplicable at large 
angles), and (b) imposing the kinematical constraint |q| < min(o;, E — u>) (in this region 
the emission angles < 1 and our formulas are qualitatively reasonable). In Fig. 2 we 
plotted udN s /du for RHIC (left) and LHC (right) conditions for (a) (upper panels) and 
(b) (lower panels) prescriptions. One can see that the energy loss is concentrated at uj < 5 
GeV. In this region the kinematical constraint suppresses the distribution by a factor ~ 5. 
It seems reasonable to assume that the results for the prescriptions (a) and (b) define an 
approximate uncertainty band of the energy loss due to the glasma color fields. 

The results presented in Fig. 2 show that the glasma should not affect strongly the 
radiative energy loss for RHIC and LHC conditions. Indeed, even without the kinematical 
constraint the distributions shown in Fig. 2 turn out to be smaller by a factor ~ 10 — 20 
than the contribution of the induced gluon emission in the QGP calculated in [9]. It is also 
seen from calculation of the total energy loss AE = J^ max du)udN s /du). Taking u m i n = m g 
and bJmax = E q /2 without kinematical constraint we obtained for RHIC (LHC) AE ^184 
(320), 234 (414), and 276 (495) MeV for E q = 10, 20, and 50 GeV, respectively. For 
the same set of E q imposing the kinematical constraint gives the values AE ^48 (63), 
49 (59), and 45 (39) MeV. These values were obtained for m g = 0.75 GeV. To estimates 
the infrared sensitivity we also performed calculations for m g = 0.4 GeV. In this case we 
obtained a rather small ~ 10 — 20% increase of AE. We also studied the sensitivity to 
the thickness L. If we choose L = 3/Q s (which is probably too large, since the glasma 
fields should become chaotic at such times) the energy loss increases by at most a factor 
~ 1.5 — 2. The obtained values of AE are clearly small compared to that from the induced 
gluon emission in the QGP which gives AE ~ 5 - 15 (10 - 30) GeV for RHIC (LHC) 
conditions at E q ~ 10 - 50 GeV [7]. 

Thus our calculations demonstrate that, despite a huge Lorentz force acting on fast 
partons in the glasma, the effect of the glasma color tubes on parton energy loss turns out 
to be rather small. This smallness is due partly to strong finite-size effects. The point is 
that, similarly to the induced gluon emission due to multiple scattering, the synchrotron- 
like gluon emission is strongly suppressed when the parton path length in the field is of 
the order of (or smaller) the gluon formation length, Lf 1 . Physically this suppression 
is connected with the fact that the initial hard parton at r = does not have a formed 
cloud of virtual gluons (with virtuality Q <E), and its splitting to a two parton state is 
simply impossible before formation of the higher Fock states in its wave function. To see 
more clearly how strong the finite-size effects are we calculated the finite-size suppression 
factor, S. We define it as a ratio S = dN s /doo j dN' s /duj, where dN s /du is our cu-spectrum 
without kinematic constraint (17), and dN s '/du is the spectrum in the approximation 
of zero formation length given by dN' s /duj = Jq dzdN s {z) / dzdou . Here dN s (z)/dzdu is 
the probability distribution of the synchrotron-like gluon emission per unit length (also 

1 The formation length of the synchrotron gluon emission for a constant f can be estimated using the 
formula L f ~ min(Li,L 2 ), where L x = 2fi/e 2 and L 2 = (24 t i/i 2 ) 1 ^ [27]. 
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without kinematic constraint) for a uniform field equals to the real local field at position 
z. In the LCPI formulation the corresponding expression reads 



dN s (z) 
dzdui 



2Re / dz 2 g[JC(p 2 ,z 2 \p 1 ,z 1 ) - JC Q (p 2 ,z 2 \p 1 ,z 1 )] , (22) 

Jzi Pl=P 2 = ' Z l= Z 



which should be calculated with the Hamiltonian (19) but for a constant f equals to the 
real Lorentz force at position z. The spectrum (22) can be written in terms of the Airy 
function [27]. For a very large slab with a smooth z-dependence of the Lorentz force, 
when at each z the local gluon formation length satisfies the inequality Lf <C z, the local 
radiation rate dN s (z)/dzdu becomes well defined. In this regime the factor S should 
become close to unity. Thus, the deviation of S from unity characterizes the strength of 
the finite-size effects. In Fig. 3 we plotted the suppression factor S for RHIC and LHC 
conditions for several quark energies. From Fig. 3 one can see that in the region u < 5 
GeV dominating energy loss the finite-size effects suppress the radiation rate by a factor 
3 — 5. Such strong suppression is not surprising since in this region the gluon formation 
length is of the order of the glasma life time. 

We presented the numerical results for gluon emission from a fast quark. For a fast 
gluon the gluon emission spectrum dN s / dxd(\ is symmetric in x -H- 1 — x. At x < 0.5 it 
is very similar in form to that for a quark, but enhanced by a factor ~ Ca/Cf = 9/4. 



4 Summary 

In this paper we studied the synchrotron-like gluon emission induced by interaction of 
fast partons with the coherent color fields in the glasma phase which should exist in 
AA-collisions at r ~ l/Q s , before formation of the QGP. We modeled the glasma by a 
time- dependent Lorentz force acting on fast partons at r < (2 — 3)/Q s . We evaluated 
the synchrotron-like gluon emission using a quasiclassical formalism based on the small 
angle approximation. But our numerical calculations show a significant contribution of 
the large angle region for gluons with u < 5 GeV which are important for the total energy 
loss. For this reason our study is not expected to give accurate information on the energy 
loss in the glasma. Nevertheless it seems relatively safe to use the results for qualitative 
estimate of the glasma effect on the radiative energy loss. 

Our results show that for RHIC and LHC conditions despite very large values of the 
coherent glasma color fields the order of magnitude of the synchrotron-like energy loss 
turns out to be small as compared to the radiative energy loss in the QGP phase. For 
this reason the glasma phase should not play a significant role in the suppression of high- 
Pt hadrons. But due to a wide angle distribution the synchrotron-like gluon emission 
may potentially be important for the jet structure at large angles (such as ridge effect and 
conical emission). However, an accurate study of this region would require the calculations 
beyond the small angle approximation. 
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Figure 1: The ratio R(q) = dN s /du)dq/ dN s /dudc{ q for q — > gq process at E q = 50 GeV 

obtained using (16) for RHIC (left) and LHC (right) conditions for the gluon energies u =2 
(solid lines), 5 (short dashed lines), 10 (long dashed lines), and 25 (dotted lines) GeV. The 
transverse momentum vector is given by q = (q, 0) (upper panels), and q = (0, q) (lower 
panels) , here x-axis in q-plane is parallel to the Lorentz force f , and y-axis is perpendicular 
to it. 
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Figure 2: The energy loss spectrum udN s /du for q — > gq process at E q = 10 (solid 
lines), 20 (short dashed lines), and 50 (long dashed lines) GeV for RHIC (left) and LHC 
(right) conditions obtained from (16) without kinematical constraint on the transverse 
momentum (upper panels) and with the restriction |q| < min(u;,.E — uj) (lower panels). 
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Figure 3: The finite-size suppression factor for q — >■ gq process for RHIC (left) and 
LHC (right) conditions for the cu-distribution obtained without kinematical constraint on 
transverse momentum, see text for details. The curves are for E q = 10 (solid lines), 20 
(short dashed lines), and 50 (long dashed lines) GeV. 
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